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I. INTRODUCTION 



We know that QCD is the right theory for the strong interactions at high energies |3J and 
also it is well accepted as the fundamental underlying theory at low energies. However, due 
to infrared slavery, it is very difficult to apply this theory in the hadronic realm because 
the resulting large coupling constant at those lower energies makes it impossible to use per- 
turbation theory. Besides ab initio calculations, like lattice techniques, one has to resort to 
an effective theory (chiral perturbation theory, ChPT ||) or sensible models which describe 
the relevant degrees of freedom restricted by the symmetries of the low energy domain. An 
example of the latter is the sigma model, || which correctly describes tttt scattering near 
threshold, to leading order in a momentum expansion (and which coincides with ChPT 
at that order). However, when one increases the center-of-mass energy, the cross section 
starts to increase until resonances are reached. There is a long standing controversy in the 
literature about the existence of an isospin zero, broad scalar resonance in tttt scattering, 
the so-called sigma meson resonance Recently, a whole conference was devoted to the 
study of the possible existence of the sigma resonance |j) . The controversy basically stems 
from the large width of the resonance, which makes it difficult to discern if the shape of 
the spectrum is actually due to a pole in the amplitude or to a result of other effects in 
the s and t channels. Several experimental results can be explained with the existence of 
such a resonance ||. If it exists, it is a relevant degree of freedom and must be incorpo- 
rated into the analysis, together with the conditions imposed by unitarity, chiral symmetry, 
etc @. It also must show up in systems other than the tttt system. For instance, it was 
pointed out ten years ago that the sigma meson can have an important role in explaining 
the AI = 1/2 enhancement in K —>■ titt decays ||. Recently it has been found strong 
experimental evidence that the sigma meson is very important in the D-meson system, in 
the singly Cabibbo-suppressed decay D + — > 7r + 7r + 7r~, being responsible for approximately 
half of the decays through the resonant sequence: D + — > tt + <j — > 7r + 7r + 7r~ ||. A best fit 
to the Dalitz plot of this decay results in m a = 483 ±31 MeV and T a = 338 ± 48 MeV, 
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where statistical and systematic errors have been added in quadrature. We want to explore 
the consequences of this experimental result in the context of the well known theories of 
light mesons and weak decays. While we were preparing this work, another study appeared 



on the effect of the a in different weak processes ||10|| . They conclude that the new result 
of E791, with a smaller value of the non-resonant contribution to T{D — > 3n) is in better 
agreement with the well measured T(D — > 2n), to which they relate via PCAC. One should 
notice however, that PCAC must be used with care, since it is valid only when the pion is 
soft, and not in the entire kinematic range of these decays. Recently, the D — > an process 
was also studied in the context of a constituent quark- meson model ||11|| . Here we study the 
specific consequences of a sigma particle and its couplings related to the D — > 3n process, 
invoking two formulations: chiral symmetry as expressed in the linear sigma model to relate 
the ann coupling to other observables and the QCD-inspired phenomenological model of 



Bauer, Stech and Wirbel [12] to treat the non-leptonic D decay. 



This paper is organized as follows. In section II we extract the effective couplings involved 
in this D decay from the experimental data. In section III we explore the consistency of 
this sigma meson within the models for spontaneous chiral symmetry breaking of strong 
interactions. In section IV we address the problem of estimating the weak process D — > an 
in the standard model of weak interactions and within the treatment of Bauer, Stech and 
Wirbel for D decays. Our conclusions are presented in section V. 



II. THE COUPLINGS IN THE DECAY D -> an -> 3vr 

From the experimental result that a fraction / = (44 ± 10)% of the D + — > n + n + n~ goes 
through the resonant channel ||, one can estimate the amplitude of the weak process D + — > 
n + a. By resonant decay in this process one understands D — > an + — > n + n~n + , i.e. a 
n + n~ pair in the final state appears through the formation of an intermediate a (sigma) 
resonance. The coupling D — a — n in this resonant process is, in general, a function of the 
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q 2 of the virtual sigma (i.e. the invariant mass of the 7r + 7r~ pair). This issue would not be a 
problem if the sigma were narrow, because then the amplitude would be strongly dominated 
at q 2 = m 2 , fixing the couplings to a constant value on the mass shell of the intermediate 
particle. As a leading approximation, we assume that the coupling is not a strongly varying 
function of q 2 and fix it on the mass shell of the sigma. A more precise estimate would 
require a model for the D decay into three pions via a sigma. Assuming constant couplings 
gDon and g a ^-n we have for the resonant three-body decay width: 



x^(l,f,f) 



1 



(q 2 —m 2 ) 2 +T a (q 2 ) 2 m 2 



where we use the standard notation A(a, b, c) = a 2 + b 2 + c 2 — 2ab — 2ac — 26c, and where 
each factor 1 / (Sn) x A 1 / 2 is the phase space integral of the corresponding two-body decay 
subprocess. The factor ^ in front is due to the symmetry of the two tt + in the final state, 
which the integral consider as distinguishable. Also T a (q 2 ) = T° x (m a /q) (p* (q 2 ) / p* (m 2 )) is 



the co-moving resonance width; here p*(q 2 ) = \jq 2 /4 — m 2 . A value for the strong coupling 
g a7Tn can be obtained from the sigma width T° by considering that the a meson decays 100% 
into Tin, two thirds of the time into charged pions. Using the general expression for the 



decay of a J=0 into two J=0 particles in terms of the coupling constant [13], we have 



2 T° a = T(a^n + 7r-)=gl 7r ^p*, (2) 
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where p* = |A 1//2 (m 2 , m 2 , m 2 ) is the magnitude of the 3-momentum of either of the final 
particles in the CM frame. 



For our numerical results we will use the data tud = 1869 MeV, m+ = 139.6 MeV [[L4J] and 
the experimentally extracted values of the a resonance (mass and width) m a = 483 ±31 
MeV and T° = (338 ± 48) MeV 0. The value for the g a7T+7T - coupling thus found is: 

g an+w . = (2.59 ± 0.19) GeV (3) 



4 



This result, together with the numerical value of the integral shown in Eq. [T], which is 
(2.55 ± 0.55) x 1CT 3 GeV~ 2 for all the appropriate values of the parameters, allows us to get 
an estimate for the weak D — a — n coupling: 

g Dan = 654 ± 120 eV. (4) 

We have thus extracted this coupling using the same prescription for the resonance that 
gave the experimental values of the sigma mass and width cited above, and also assuming 
that the coupling is independent of the virtuality of the intermediate a state. The validity 
of our estimate is thus conditioned to the assumption that the coupling for q 2 values away 
from the resonant peak does not differ by much from its value at the peak and that the 
shape of the resonance away from the peak is as prescribed in Eq. [TJ Both approximations 
are clearly related and constitute the old issue of the background underlying the resonance. 

To make a more crude estimate that does not deal with the background issue, one could 
use a narrow width approximation to estimate the coupling. Although this is a priori not a 
very good approximation since the a resonance is not narrow, the estimate is rather robust 
because it hides the width of the resonance and only deals explicitly with its branching ratio. 
Within this approximation the resonant D decay width is: 

T(D + -> ir + ir + iT-) res « \ x T(D + -> a7r + ) x Br(a -> tt+tt"), (5) 

where again there is a symmetry factor 1/2 due to the two identical 7r + . Now, Br [a — > 
7r + 7r~) = 2/3 and T(D + — > 7i + n + 7i^) res = / X T(D + — > 7r + 7r + 7r~), according to the experi- 
mental analysis. Using the established data Br(D + — ► 7r + 7r~7r + ) = (3.6 ± 0.4) x 10~ 3 and 
t d + = (1.06 ± 0.02) x 10~ 12 s [fL4H , one thus obtains the estimate for the decay D + — ► an + in 
this approximation: 

Y(D + -»■ a7r + ) = 3 x / x Y(D + -> = (2.94 ± 0.75) x 10" 12 MeV. (6) 



From this rate we can again extract the gnaw coupling using the general expression for the 
decay of a J=0 into two J=0 particles: 
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(7) 



where p* = l/2x A 1 / 2 (mf^, m 2 , m 2 ) is again the 3-momentum of either of the final particles 
in the CM frame. We thus obtain 



One can notice that the central value of this estimate is only 16% smaller than our more 
elaborate estimate of Eq. |j. In this sense, however crude the narrow width approximation 
could be considered a priori, it is indeed quite robust, since the result for gnan thus obtained 
does not differ much from our more elaborate estimate, in which the full resonance shape is 
taken into account. 



Having estimated the couplings directly from the experimental data, let us see the conse- 
quences of these results within the theory. First, consider the aim coupling. We are here in 
the low energy regime of strong interactions. Consider then the linear sigma model as the 
framework, with the a meson as the scalar (J PC = ++ ) that remains from the breakdown 
of chiral symmetry. It is sufficient to consider QCD with the lightest quarks u and d only, so 
that the chiral symmetry in question is SU (2) l x SU(2)r broken down to SU (2)y or isospin. 
The meson sector comprises the massive a scalar and the three pions (pseudoscalars) 7r a , 
a = 1, 2, 3, which play the role of Goldstone bosons. Using the notation of ref. |l5j the quar- 
tic coupling of the potential, A, and the vacuum expectation value of the field, v, determine 
all the physical parameters, like the mass of the sigma as m a = 2Xv 2 , the aim coupling as 
g ann = 2Xv and the (charged) pion decay constant as f n = \/2v. The linear sigma model 
thus predicts the g U7T7T coupling in terms of the mass of the a particle and the pion decay 
constant: 



gocj-K 



548 ± 70 eV. 



(8) 



III. THE a MESON IN THE LINEAR SIGMA MODEL 
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g awn = ^ = (2.54 ± 0.01) GeV, (9) 

which is surprisingly close to the value deduced from the data on resonant D decay. There- 
fore, there is an apparent consistency of the linear sigma model description for the 
resonance in D decays. However, a word of caution is due here: the validity of a pertur- 
bative calculation within a linear sigma model treated as a quantum theory is questionable 
if the mass of the scalar o is large (i.e. if m a ~ /„.), which is precisely the case in the 
real world. A perturbative calculation could be trusted if the expansion parameter is much 
smaller than unity. In the sigma model, the expansion parameter is A, or more precisely 
6A/(167r 2 ). From the data on m a and f n it is easy to see that the value for the expansion 
parameter is: 

& ~ 0.5, (10) 

which is uncomfortably close to unity. This is precisely the reason why one would prefer 
to use a non-linear version of the model as a perturbative theory, where the heavy sigma 
is integrated out and only the light particles (the pions) are considered. The effective 
interaction among the pions is then explicitly weak at low momenta, vanishing at threshold 
(a feature which is evidence that chiral symmetry is spontaneously broken at low energies 
and that the pions are the corresponding Goldstone bosons). Here we used the linear sigma 
model with the sole purpose of exhibiting the sigma particle in the theory. 

IV. THE D+ -> CJ7T+ FORM FACTOR IN THE BSW MODEL 

The coupling goan found in Eq. |], which has dimensions of mass, is an effective result due 
to an underlying weak interaction process. In this section we want to examine this coupling 
at that fundamental level. We therefore estimate the relevant matrix elements of weak 
interactions that contribute the D + — > <jtt + decay. However, because of the effect of strong 
interactions, the asymptotic states are not the elementary quarks but composite hadrons, 
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meaning that the analysis cannot be done purely at the level of fundamental interactions 
treated perturbatively. Only the weak interaction can be considered at leading order, but 
to take into account strong interactions we require of a sensible model. We thus study the 
problem within the well accepted model of Bauer, Stech and Wirbel (BSW) for non-leptonic 
D decays |T2|] . For our purposes, the a particle will be a J PC = ++ isoscalar with quark 
content (uu+dd). In order to tackle non-leptonic decays, one has to use the operator product 
expansion to construct an effective weak hamiltonian containing local four-quark operators 
with Wilson coefficients that can be computed perturbatively at the appropriate energy 
scale fi. In our case, the relevant effective hamiltonian is: 

H w = ^V* d V ud {cx(li){dc)(y- A ){ud)(y- A ) + C 2 (/j)(uc)(y-A)(dd)(y-A)} , (H) 

where (qq')(v-A) = <tin (1 _ 7s) <?'• In the BSW approach,the effective weak hamiltonian 
density is built in terms of products of hadron currents which mimic the underlying quark 
currents of the weak interactions: 

H eff = %V: d V ud { ai J^J^ + a 2 J^J^} . (12) 

The coefficients ai and a 2 that accompany the products of currents are now phenomenological 
constants of the model that must be fitted from experiment. In this approach, the decay 
amplitude factorizes into a product of two current matrix elements. We will use this effective 
hamiltonian to study the D + — > an + decay process. 

The factor proportional to a x in the matrix element (o"K + \H e ff\D + ) involves two terms. 
One of them is (Vtt + | |0^} (o\ J A dc ^ |-D + ^>, so-called "annihilation" term, which we neglect 
within the model as it is prescribed, because it involves form factors of the light mesons 
at high momentum (q 2 = m 2 D ) fl2|| . We should emphasize here that we are not neglecting 
the so-called quark annihilation diagrams that play a role in explaining D meson lifetime 
differences and other controversial issues, but we are only following a consistent prescription 
of the model. Indeed, in the BSW model the coefficients a\ and a 2 are only phenomenological 
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parameters and do not exactly correspond to the Wilson coefficients ci(/i) and C2(/i) obtained 
in perturbative QCD and associated with the quark amplitudes. The other term, which is 
the only relevant term here, corresponds to (cr(p')\ J^ dc ^\D + (p)^ (n + (q)\JA^\®)- After this 
factorization, one needs the matrix element of the current (notice that only the axial currents 
J a contribute, because of the parity of the hadrons involved): 

(<T(p')\jf c) \D+(p)) = F^\q 2 ) [{p + p'T - ^^/) + F^\q 2 )^^lq^ (13) 

with -Fi(O) = Fq(0) in order to avoid a spurious singularity at q 2 = 0. The two form factors 
Fi and Fq correspond to transverse and longitudinal components of the current, respectively. 
If these form factors are dominated by poles, those of Fi should be at axial vector meson 
masses and those of F at pseudoscalar meson masses, all of cd flavor content. However, 
only F enters in our case, because (7r + (g)| J^|0) = if-^cf is purely longitudinal. 

The factor of 02 in the matrix element (a7i + \H e ff\D + ) vanishes, because it involves (a\ J M |0) 
which is identically zero due to conservation of the vector current. Therefore, from all of 
the above, the D-a-n coupling is: 

an+\H eff \D + ) = g Da7T = ^V^V ud a, F^ (m 2 ) x (m 2 D - ml) /„. (14) 



Using V* d V u d = 0.21, the model value a± = 1.10 ± 0.05 fitted for D decays [16| and our value 
for the coupling gDan i n Eq. ^, we get an estimate for the axial form factor 

F o Da \™l) = 0.79 ±0.15 . (15) 

We can safely extrapolate from q 2 = m\ to q 2 = 0, since we do not expect poles at light 
masses for a charmed current: 

F o (0) = Fi(0) w 0.8 ±0.2 , (16) 



which agrees within errors with the form factor of Ref. |T3] for D — > tt and is within the 
range of most of the other D decays (0.6 to 0.8). A remarkable fact of this result is that 
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the form factors in Ref. [12| are calculated assuming the mesons are bound states of the 



corresponding valence quarks, while in our case one would hardly treat a priori such a short 
living resonance like the sigma as a bound state, let alone using a wave function for it. As 
expected, strong phases do not play a role in the D + — ► o~7r + decay since the final state has 
only one definite value of isospin ,1—1. Strong phases in the (3ir) final state have already 
been taken into account by the introduction of resonances in the fitting procedure to the 
Dalitz plot. 

V. CONCLUSION 

We have used the new experimental evidence for the a meson in the non-leptonic D + — > 
7r + 7r + 7r~ decay to estimate the effective goa-w and g an7T couplings from the data and the 
resonance shape. We then studied the consequences of the values of these couplings within 
the accepted theoretical framework and found consistency of the latter with the experimental 
data. 

Acknowledgments 

We would like to thank Ignacio Bediaga for the suggestion of and useful comments on this 
work and to Jon Rosner for reading the manuscript. One of us (R.R.) would like to thank 
the hospitality of the Departamento de Fisica of the Universidad Tecnica Federico Santa 
Maria, Valparaiso, Chile, where this work was carried out. This work was supported in part 
by CNPq, Pronex (Brazil) and FAPESP (Sao Paulo) and by Fondecyt, Chile, grant No. 
1980150. 



[1] See, e.g., I. Hinchliffe and A. V. Manohar, |iep-ph/ 0004186 for a review. 



10 



[2] See J. Gasser, Nucl. Phys. B279, 65 (1987) and references therein. 

[3] M. Gell-Mann and M. Levy, Nuovo Cimento 16, 705 (1960); S. Weinberg, Phys. Rev. Lett. 18, 
188 (1967). 

[4] V. Elias and M. Scadron, Phys. Rev. Lett. 53, 1129 (1984); R. Kaminski, L. Lesniak and J.P. 
Maillet, Phys. Rev. D50, 3145 (1994); N. T6rnqvist, Z. Phys. C68, 647 (1995); N. T6rnqvist 
and M. Roos, Phys. Rev. Lett. 76, 1575 (1996); ibid Phys. Rev. Lett. 77, 2333 (1996) and 
Phys. Rev. Lett. 78, 1604 (1997); N. Isgur and J. Speth, Phys. Rev. Lett. 77, 2332 (1996); M. 
Harada et al, Phys. Rev. D54, 1991 (1996); ibid Phys. Rev. Lett. 78, 1603 (1997). For more 



recent views on the problem, see M.R. Pennington, preprint tiep-ph/9905241, talk given at 
Workshop on Hadron Spectroscopy (WHS 99), Rome, Italy, 8-12 Mar 1999; S. Ishida, preprint 
[hep-ph/9905260 , talk given at Workshop on Hadron Spectroscopy (WHS 99), Rome, Italy, 



8-12 Mar 1999; L. Celenza, H. Huang and C. Shakin, CUNY preprint BBCCNT 00/051/290 
(May 2000), submitted to Phys. Rev. C. 

[5] See, e.g., N.A. Tornqvist, summary talk at Conference on the Possible Existence of the Light a 
Resonance and its Implications to Hadron Physics, Yukawa Institute for Theoretical Physics, 
Kyoto, Japan, ll-14th June 2000, [hep-ph/0008135| . 



[6] J.E. Augustin et al. (DM2 coll.), Nucl. Phys. B320, 1 (1989); C. Amsler et al, (Crystal Barrel 
coll.), Phys. Lett. B291, 347 (1992); ibid B322, 431 (1994); B342, 433 (1995); B353, 571 
(1995); F. Antinori et al. (WA91 coll.), Phys. Lett. B353, 589 (1995); S. Ishida et al, Progr. 
Theor. Phys. 95, 745 (1995); M. Ishida et al, Progr. Theor. Phys. 101, 601 (1999); ibid, 
preprint Eiep-ph/9905261 , to appear in proc. of Workshop on Hadron Spectroscopy (WHS99), 



Rome, Italy, 8-12 Mar 1999. 

[7] E. van Beveren et al, Z. Phys. C30, 615 (1986); J. Oiler et al, Nucl. Phys. A620, 438 (1997); 
ibid Phys. Rev. D59, 074001 (1999) and D60, 074023 (1999); R. Kaminski et al, Eur. Phys. 
J. C9, 141 (1999); M. Scadron, Eur. Phys. J. C6, 141 (1999). 



11 



[8] T. Morozumi, C.S. Lim and A.I. Sanda, Phys. Rev. Lett. 65, 404 (1990). 

[9] J. Miranda, talk presented at the XX Meeting of the Brazilian Physical Society- Session Parti- 
cles and Fields, October 25 - 28, 1999; A.J. Slaughter, talk presented at the XXXVth Rencon- 



tres deMoriond, March 18 - 25, 2000; E791 Collaboration, E.M. Aitala et al, [hep-ex/0007028 
submitted to Physical Review Letters. 



[10] A.D. Polosa, N. A. Tornqvist, M.D. Scadron and V. Elias, preprint |hep-ph/0005106 . 

[11] R. Gatto, G. Nardulli, A.D. Polosa and N. A. Tornqvist, Phys. Lett. B494, 168 (2000). 

[12] M. Bauer, B. Stech and M. Wirbel, Z. Phys. C34, 103 (1987). 

[13] See, e.g., J.L. Rosner, Phys. Rev. D60, 114026 (1999). 

[14] C. Caso et al. (Particle Data Group), Eur. Phys. J. C 3, 1 (1998). 

[15] M.E. Peskin and D.V. Schroeder, An Introduction to Quantum Field Theory (Addison- Wesley 
Publ. Co., 1995), page 349. 

[16] M. Neubert and B. Stech, in Heavy Flavours II, edited by A.J. Buras and M. Lindner (World 
Scientific, 1998), page 294. 



12 



